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1 Introduction 

Lattice simulations in the Landau or maximal Abelian gauges have made clear the 
need for mass parameters in the fitting functions [H [21 [3] . These fits involves e.g. a 
Yukawa propagator ~ g2_|^j^^2 • Moreover, phenomenological studies frequently make 
use of an effective gluon mass [1]. Theoretical studies, based on Schwinger-Dyson 
equations and the Pinch Technique report a dynamical gluon mass [3, El [7]. Also 
alternative analytical calculations gave such evidence, see e.g. [8]. Finally, let us 
mention the issue of ^ power corrections in physical correlators, tackled with QCD 
sum rules lattice and OPE techniques [10] or even via the AdS/QCD picture 
[TT] . A natural question arising is where does such a mass scale originate from? We 
recall that the standard Yang-Mills action cannot contain a mass due to gauge in- 
variance. The Higgs mechanism is not a solution here, due to the associated gauge 
symmetry breaking. A natural answer is that a dynamical mass scale is generated by 
nonperturbative effects, in the form of a dimension 2 condensate. Of course, then the 
question pops up which d = 2 operator O to consider? In our opinion, a few require- 
ments are to be met: O should be gauge invariant, as it is supposed to enter physical 
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quantities. It should be local, as nonlocal actions are hard to handle/interpret, and 
it should be renormalizable, as we want to perform quantum calculations, thus we 
want consistent finite results, renormalization group (RG) invariance, .... Zakharov 

but nonlocal [12]. Only in the Landau gauge do we have (^min) = ^ which is a 
renormalizable and condensing local composite operator, giving rise to effective dy- 
namical gluon mass [H [13] . We could also use the mass operator based on the trans- 



verse gluon field: (A^)^ 



U = 



A,. 



n 2 



or the non-Abelian Stueckelberg term. 



The previous proposals are all classically 
but face problems with nonlocality, nonpolynomiality, nonrenormaliz- 

ab 



Os = {a, 
equivalent 

ability, .... We proposed to study O 
to the Euclidean YM action via 

_ 



arpa 



F^u ' which we couple 



ab 



where D^r denotes the covariant derivative. This operator O already found use in 3D 
YM p5]. This operator is clearly gauge invariant. Furthermore, it can be put quite 
easily in a local form, since we can replace So with 

S'o - \ j d'^Bl^Dl'Dl^Bl^ + \ J d'xGl^Dl'D'^G'^^^ + '^jd'^^B- 3)^ F;, , (2) 

where B'^^, B^^^ are antisymmetric bosonic and G"^^,G°'^^ fermionic (ghost) fields in 
the adjoint representation. Notice that for m = 0, we have in fact introduced a unity. 



2 Renormalization analysis 

We continue our investigation from the starting action 

V ' V ' 

SyM Sgf 

+ I J d^xGlD^'Dl^G^^^ + lI^Jd^,^B- B)l f;, . (3) 

As the reader will notice, we imposed a linear gauge fixing, encoded in Sgf. The 
complete action enjoys a nilpotent BRSTi symmetry, generated by 

s^A1 = -Dfc', s,^^^-r''c'c\ s^B^, ^ g B^, , s^B;^ ^ gr'^c'B;^ , 

siGl, = gr'^c^G^. , siGl^^gr'^c'Gl,, s{e s.b^ ^ , = . (4) 

For m = 0, we can identify a nilpotent "supersymmetry" Sg |16] . 

(5,5^, = G;:, , 6sK. = B", 5,^ = for all other fields ^ , 5^ = 0, (5) 
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which can be invoked to estabhsh another nilpotent BRST2 invariance for m = 0, 



S2 = Si + 5s 



0. 



(6) 



For m 7^ 0, we can embed the action ^ into a "larger" model with external sources 
[14j . In a particular physical limit, the original model is recovered. The reason 
why we introduced the "larger" model is because it simplifies the renormalizability 
analysis. If we can prove the renormalizability of this action, then also the physically 
relevant action will be renormalizable as a special case. The extended model obeys 
many Ward identities, including a Slavnov- Taylor identity. We constructed the most 
general action compatible with these Ward identities, and when the smoke cleared 
after taking the physical limit, we arrived at a renormalizable local action with mass 
terms [T6l. 



Sphys 
Sd 



Scl + Sgf , 
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a 



+ 



^abcd 



B 



16 



^ p<T^ pa 



(7) 



which is still BRSTi invariant (but not BRST2\), whereby the classical part Sd is 
gauge invariant. For reasons of renormalizability, we had to introduce the scalar 
(mass) couplings Ai^s and the gauge invariant tensor coupling X"-'"^'^, Notice that the 
new quartic interaction oc A'^*^'' in the novel fields spoil the unity. One might therefore 
question the equivalence with massless YM theories when m = 0. In [16], we have 
been able to show that 



(YM functional) 



Sym+S, 



gf 



(YM functional) 



Sphys+Sgf 



by making use of the ^^-cohomology, meaning that expectation values of "pure" YM 
functionals remain unchanged when calculated in our massless model. As a conse- 
quence, there will be no X"-^'^'^ independence in those pure YM Green functions, and 
RG functions of the original YM quantities will remain unchanged, as long as mass- 
less renormalization schemes are employed. This fact was confirmed by explicit 1- 
and 2- loop computations in [TU [16]. In addition, also several other RG functions 
were computed. The obtained results were consistent with the Ward identities and 
confirmed the renormalizability explicitly. We end this section by mentioning that 
although the BRST2 symmetry is broken, with an associated broken Slavnov- Taylor 
identity ST2, it is nevertheless possible to show that the ensuing Ward identities 
between a large class of Green functions are as if the ST2 would be unbroken [T7j . 
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3 Unitarity analysis 



Since we have a massive gauge model with nilpotent BRSTi symmetry generator, we 
might hope that the theory would be unitary. Assuming that we start from the action, 
and that we take the elementary field excitations as asymptotic states, how can we 
prove the unitarity of the iS-matrix? This rather complicated task was performed in 
|18j . whereto we refer for all details. One of the main tools was the use of the (free) 
BRST cohomology. We start from the free action 5*0 with free BRST symmetry sq 
with nilpotent charge Qo {Ql = 0). Physical states It/jp) are defined as belonging to 
the Qo-cohomology, i.e. 

I^p) e nj,Hys ^ Qol^p) = , \ijp)^\...) + Qo\...), Qol • ■ .) = . (9) 

Then the 2 remaining questions are: is this definition invariant under time evolution, 
described by the 5-matrix, and do the physical states have a positive norm, a conditio 
sine qua non for a sensible quantum theory. Concerning the time evolution, we recall 



that in the operator language S = T 



so that Sl'ipp) G Hphys, if we 



require that [S, Qo] = 0. We can pass to the path integral language to rephrase this 
condition into one for the action S, namely 

sqc'"^ = on-shell, i.e. modulo free equations of motion . (10) 

This equation can then be solved iteratively, in particular order by order in the 
available coupling constant (s). 

action S = Sq + Si + . . . , BRST s = Sq + Si + . . . , 

while it is proven that 

{so + ... + Si){So + ... + S^) = 0, (so+ ... + SiY = to i"" Older. (11) 

If we are lucky, the procedure stops at finite order, and we end up with an action 
S, invariant under a BRST symmetry s with nilpotent generator Q, with the desired 
property that the physical subspace Tiphys is invariant under time evolution. We used 
a "backward" argument: if we start from an action S with nilpotent symmetry charge 
Q, then this is a solution of the previous procedure starting from the corresponding 
free counterparts .So and Qo, obtained by switching of any couplings. Having answered 
the time evolution question, we are still left with the positive norm issue. We suffice 
by saying this required a rather lengthy and technical study. One of the problems 
faced was the occurrence of multipole fields. Nevertheless, a complete analysis was 
provided in [T8], yielding a negative result: negative norm states do appear in the 
physical subsector of the theory, therefore it is not unitary. The reader might have 
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immediately questioned the wisdom of trying to prove the unitarity of this model, 
as it has its root in a nonlocal field theory, which are known to have problems with 
ghost states. Next to the presence of a nilpotent BRST invariance, there is however 
another reason why this endeavour was worth the effort. Specializing to the Abelian 
case, it can be shown that after integrating out the auxiliary fields, we obtain the 
Abelian Stueckelberg model, where the Stueckelberg scalar has been integrated out. 
We recall that the Abelian Stueckelberg model is renormalizable and unitary, see e.g. 
[T9] for a review. If we would analyze the unitarity of the Abelian version of our 
model, we would run into exactly the same problem as in the non-Abelian case, i.e. 
the presence of negative norm states in the physical subspace. Apparently, the way 
of localizing a nonlocal action plays a substantial role. 

4 Discussion 

Since our action is perturbatively equivalent with YM in the massless case, we could 
take it as starting point instead of YM. As proven, we can couple mass terms to it, 
without ruining renormalization/gauge invariance requirements. The next challenge 
would of course be trying to construct a sensible gap equation to produce a dynami- 
cally generated value for m oc Aqcd- Afterwards, we could start looking at potential 

2 

nonperturbative ^ power corrections appearing in gauge (in)variant correlators. We 
are thus interested in nonperturbative effects in an asymptotically free theory, which 
occur in an energy region below the high energy (asymptotic) region where the el- 
ementary fields excitations are observables. Since this is YM (m = 0), there is no 
unitarity issue at high energies: 2 transverse gluon polarizations are physical. At 
lower energies, the gluons etc still are our effective degrees of freedom, but will be- 
have like quasi particles, corrected by nonperturbative effects due to the generated 
m 7^ 0. Lack of unitarity in terms of gluons is hence not a problem, but rather finding 
the correct physical degrees of freedom, which of course corresponds to the task of 
proving confinement. 

Acknowledgments. 

D. Dudal is a Postdoctoral Fellow and N. Vandersickel a PhD Fellow of the Research 
Foundation-Flanders (FWO Vlaanderen). The Conselho Nacional de Desenvolvi- 
mento Cientifico e Tecnologico (CNPq-Brazil), the SR2-UERJ and the Coordenagao 
de Aperfeigoamento de Pessoal de Nfvel Superior (CAPES) are gratefully acknowl- 
edged for financial support. 



5 



References 

[1] K. Langfeld, H. Reinhardt and J. Gattnar, Nucl. Phys. B 621, 131 (2002). 

[2] K. Amemiya and H. Suganuma, Phys. Rev. D 60, 114509 (1999). 

[3] V. G. Bornyakov, M. N. Chernodub, F. V. Gubarev, S. M. Morozov and M. I. Po- 
likarpov, Phys. Lett. B 559, 214 (2003). 

[4] J. H. Field, Phys. Rev. D 66, 013013 (2002). 

[5] J. M. Cornwall, Phys. Rev. D 26, 1453 (1982). 

[6] A. C. Aguilar and A. A. Natale, JHEP 0408, 057 (2004). 

[7] A. C. Aguilar and J. Papavassihou, JHEP 0612, 012 (2006). 

[8] H. Verschelde, K. Knecht, K. Van Acoleyen and M. Vanderkelen, Phys. Lett. B 
516, 307 (2001). 

[9] K. G. Chetyrkin, S. Narison and V. L Zakharov, Nucl. Phys. B 550, 353 (1999). 

[10] Ph. Boucaud, A. Le Yaouanc, J. P. Lcroy, J. Micheli, O. Pene and J. Rodriguez- 
Quintero, Phys. Rev. D 63, 114003 (2001). 

[11] O. Andreev, Phys. Rev. D 73, 107901 (2006). 

[12] F. V. Gubarev and V. L Zakharov, Phys. Lett. B 501, 28 (2001). 

[13] D. Dudal, H. Verschelde and S. P. Sorella, Phys. Lett. B 555, 126 (2003). 

[14] M. A. L. Capri, D. Dudal, J. A. Graccy, V. E. R. Lcmcs, R. F. Sobreiro, 
S. P. Sorella and H. Verschelde, Phys. Rev. D 72, 105016 (2005). 

[15] R. Jackiw and S. Y. Pi, Phys. Lett. B 403, 297 (1997). 

[16] M. A. L. Capri, D. Dudal, J. A. Gracey, V. E. R. Lemes, R. F. Sobreiro, 
S. P. Sorella and H. Verschelde, Phys. Rev. D 74, 045008 (2006). 

[17] M. A. L. Capri, D. Dudal, V. E. R. Lemes, R. F. Sobreiro, S. P. Sorella, R. Thibes 
and H. Verschelde, Eur. Phys. J. C 52, 459 (2007). 

[18] D. Dudal, N. Vandersickel and H. Verschelde, Phys. Rev. D 76, 025006 (2007). 

[19] H. Ruegg and M. Ruiz-Altaba, Int. J. Mod. Phys. A 19, 3265 (2004). 



6 



